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A simple model is put forward which accounts for the 
occurrence of certain generic dewetting morphologies in thin 
liquid coatings. It demonstrates that by taking into account 
the elastic properties of the coating, a morphological phase 
diagram may be derived which describes the observed struc- 
tures of dewetting fronts. It is demonstrated that dewetting 
morphologies may also serve to determine nanoscale rheolog- 
ical properties of liquids. 



When a thin hquid fihn beads off a solid substrate, 
it is eventually transformed into an ensemble of individ- 
ual liquid droplets, the arrangement of which may vary 
strongly according to the basic mechanisms involved in 
the dewetting process In by far most cases, dewet- 

ting is initiated by heterogeneous nucleation of individ- 
ual holes in the initial film, thus forming contact lines 
between the film surface and the substrate. The surface 
forces acting upon these contact lines give rise to mov- 
ing dewetting fronts, the motion, shape, and interplay of 
which largely determine the final dewetting morphology. 

It would be very desirable to have a theory describ- 
ing the dynamics of these fronts precisely enough for ex- 
tracting information on the dewetting mechanisms from 
the characteristics of their shape and the final droplet 
structure. However, the task of setting up such a theory 
has proved untractable so far. Solving the Navier-Stokes 
equation with the moving contact line as a boundary con- 
dition is particularly complicated and entails ad hoc as- 
sumptions on the dynamics at the contact line [p|-p^. 
Furthermore, it has meanwhile become clear that vis- 
coelastic effects, which are not described by the Navier- 
Stokes equation, may be decisive for the evolving mor- 
phology in both advancing and receding fronts [ lO-Q . 
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In particuilar, it has recently been shown that 
'elastic dewetting fronts may strongly differ, as to their 
shape and dynamics, from Newtonian ones ||l5[] , and it 
has been pointed out that these findings are not ac- 
counted for by theory so far Consequently, what 
is needed is a handy theoretical model, which is just as 
detailed as necessary to describe the main physics in- 
volved in the dewetting dynamics of viscoelastic films. It 
is the purpose of this paper to propose such a model. 

Let us first recall the main features to be explained, by 
summarizing the generic front morphologies observed so 
far. Fig. ^ shows dewetting fronts in liquid polystyrene 
films beading off silicon substrates (the fronts move from 



left to right). The profiles are obtained by scanning force 
microscopy from the rims of circular dry patches nucle- 
ated in the films |p6|-p^. Most commonly, one observes 
profiles as those shown in fig. with a simple decay 
from the crest into the undisturbed film to the right. 
However, when the molecular weight of the polymer is 
small enough , the shape may be qualitatively differ- 
ent. This is shown in fig. which has been obtained 
with polystyrene films with a molecular weight of only 2 
kg/mol. A ditch is clearly visible in front of the crest, and 
even another small elevation to the right of the ditch is 
present as well, such that the front appears as a damped 
oscillation. The only difference between the film materi- 
als used for fig. [|a and |l|b is their molecular weight, and 
thus their viscoelastic properties. If the film is not too 
thick, the ditch may reach the substrate and pinch off 
the crest from the film, forming a new contact line. This 
happens repeatedly, such that a series of isolated crests 
forms in a cascade of pinch-off events, as shown in fig. |l|c. 

We shall now create a suitable mathematical descrip- 
tion of these findings. An incompressible viscoelastic 
fluid may be described by the force balance Jl9|-pl| 



(1) 



where rj is the viscosity, j is the material current in the 
film, p is the pressure, E is Young's modulus, and 4> is 
a local displacement field describing the strain pl| , p2t . 
In the present treatment, we will adopt the lubrication 
approximation ]2^ ], thereby neglecting the normal {z-) 
components of the current and the pressure gradient, and 
consider quantities averaged over the film thickness, h. 
We then can rewrite eq. (|l]) as J{x) — —C'dx{p — aEdx(f>), 
where J is the total current in the film, and a is a nu- 
merical factor which characterizes the flow profile. The 
latter depends, e.g., on the friction of the film at the sub- 
strate, and it is clear that in the case of full slippage, the 
coupling of the flow to the strain will not be the same 
as in the absence of slip. These effects are difficult to 
treat explicitly, and thus will be absorbed here into a. 
Similarly, C depends upon the viscosity of the film, 77, 
and its friction coefficient at the substrate surface, k. If 
the latter is infinite (Poiseuille flow with no slip at the 
substrate), we have C — h? /irj, while for k ^ (plug 

flow), c= /iVk mn]. 

Since the derivative of the displacement is the strain, 
a{x,t) = 3x4', s-nd excess pressure in the film comes 
about from the curvature of its surface, we have 



J{x, t) = Cdx{ldxxC,{x, t) + aEa(x, t)} 



(2) 
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where 7 is the surface tension, and ({x,t) is the ver- 
tical displacement of the film surface. We explicitely 
neglect long range wetting forces acting through the 
film, such as van der Waals forces. In sufficiently thick 



films, which dewet by nucleation |17 24[| (such as assumed 
above), these do not play a significant role. 

Since the strain in the film may decay by internal re- 
laxations of the material according to dtcr — —luoct, its 
coupling to the material current can be written as [pl|j22| 



{dt + ujQ)a = dxJ/h 



(3) 



where we have neglected the nonlinear convective term 
Jdx4>- This restricts our discussion to small excursions 
C,, but enables a linear treatment. By combination of 
the above equations with the continuity equation, dtC + 
dxJ = 0, it is easy to obtain the equation of motion for a 
and ^. A precondition for elasticity to play a significant 
role is that wq << cj, in which case the dispersion relation 
reads 



^ 2 / 2 
OJ = LOQ + Cq \ — 



(4) 



where q is the wave number of the perturbation. 

In order to determine the shape of the moving dewet- 
ting front, we now look for travelling- wave solutions of 
the form C,{x, t) — C,{x — vt), where v is the velocity of the 
dewetting front. We remain in Fourier space, i.e., we de- 
compose the front profile into modes Qq oc expiiqx — cut). 
Writing Q :~ iqh, it follows that uj — Qv/h, and from 
eq. (||) we then get 



aEh 



7 



Q + — — 



0. 



(5) 



As it will become clear below, the term cooh'^/jC is rather 
small. If we neglect it for the moment, we get the cubic 
equation 



(6) 



where we have introduced the dimensionless quantities 



V := ^ and E := This can be easily solved, and 

the solutions of the quartic equation (H) turn out to devi- 
ate significantly from those of the cubic equation only for 
parameters which are not relevant in most experimental 
situations. We will thus primarily discuss eq. (^, and 
refer to eq. only occasionally, where appropriate. 

If viscoelastic effects are absent, i.e. E = 0, we have 
Q = v^. We are interested only in modes with ^{Q) < 
p5[ , which are Q — — iv^(l ± The shape of the 

real surface of the leading front is thus oc cos{^/3x) 
0,|6|, which clearly exhibits not only a ditch, but a 
damped oscillation. For Newtonian liquids, it can be 
seen from numerical solutions of the Navier Stokes equa- 
tion for the moving front problem that this result does 



not depend on the friction of the film at the substrate 
[|ll|,|2|. A profile like that is shown in fig. ||b, and it is 
obtained only for small molecular weight indeed, where 
elastic effects are particularly small. The real and imag- 
inary part of Q can be determined from the profile, and 
in the case of the one displayed in fig. it turns out 
that they are not precisely as expected for a Newtonian 
fluid iQ . Since the complex roots of eq. (|^) satisfy the 
relation 



35R(Q)2 ^ 3(Q)2 = E 



(7) 



the elasticity of the film can be inferred from the mea- 
surements, and we get aE = 8.7 kPa. We would like 
to point out that this is a non-invasive method of deter- 
mining rheological properties of liquids on a sub-micron 
scale. 

As the height of the crest, H, increases while more 
material is collected in the front, the amplitude of the 
oscillation increases, and so does the depth of the ditch 
Q. When the latter reaches the substrate, the crest is 
pinched off and a cascade is formed as shown in fig. |l|c. 
However, this works only as long as the width, W, of 
the crest is less than half the wavelength of the damped 
oscillation. As H increases, so does W |17 2^, with the 
ratio H/W =: G(0) being only a function of the dynamic 
contact angle at the substrate, 8 jl^. For a cylindrical 
rim [|3|, we have simply G{Q) = (1 — cos6)/2sin6. 
For more asymmetric shapes, the form of G is different, 
but similar as to the overall behavior and the order of 
magnitude. If W is larger than Wc = A/2 = 7r/3(Q), the 
depth of the ditch is rather determined by the contact 
angle of the rim at the leading edge. Consequently, it 
will reach a final maximal depth iQ, and a dry spot 
forms only if this depth exceeds the film thickness. The 
condition for a cascading front to form is readily seen to 
be 



7rG(e) >| 3(Q) I -exp 



(8) 



As the film thickness is decreased, so is the right hand 
side of eq. (||) {h cancels out in the exponential) , such that 
a cascade is expected at sufficiently small h. Eq. (^) can 
be used to determine the boundary line for the appear- 
ance of cascade structures in the {v, ii^)-plane. In fig. ||a, 
this is indicated by the lower curve for G = tt^^. For 
different values of G, the maximum of the boundary line, 
indicated by the dot in the figure, is quite accurately de- 
scribed hy E = 0.060 • G^. For very thin films, where 
van der Waals forces play a role, the cascade region will 
slightly extend or shrink, depending on the sign of the 
Hamaker constant. 

Let us now turn to the infiuence of stronger elastic 
effects. It is interesting to investigate at which system 
parameters the damped oscillation, as shown in fig. 
vanishes. These are given by the zero of the discriminant 
of the cubic equation (^), which is 
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(9) 



The critical modulus above which there is no oscillation 
is thus given byi? = 3(|)^/^ This may be viewed as an- 
other morphological 'phase boundary' in the(i?,{;)-plane, 
and is shown in fig. ^ja as the upper solid curve. The exact 
boundary, according to eq. (|^), deviates noticeably only 
for very small values of E and v. To quantify this devia- 
tion, it is useful to introduce the dimensionless parameter 
n := LOoh/v, which represents the importance of intrin- 
sic relaxation in the material with respect to the shear 
flow exerted by the dewetting process itself. In fig. ^d, 
the 'exact' phase boundary, which was determied by solv- 
ing eq. (H) graphically, is shown as the solid curve, along 
with the approximation obtained from the cubic equation 
(dashed line). For 17 — > 0, the solid curve approaches the 
dashed line asymptotically. The diagram presented in 
fig. which corresponds to = 0, shall thus be consid- 
ered quite reliable within the range displayed if < 0.1 
(in the experiments in ref. Q w 10~^). 

The two complex conjugate roots giving rise to the os- 
cillations for D > merge into a degenerate real root 
at a critical elastic modulus {D — 0). As the latter in- 
creases further, they bifurcate again into two real roots, 
Qiong and Qshort- In fact, 'elastic' fronts can be well fit- 
ted by a superposition of two exponentials, as shown in 
fig. |l|a by the dotted curves. The model predicts the ratio 
Q short IQiong to incrcasc with and thus the front to 
become more and more asymmetric, in accordance with 
observation 13 1^. We can again determine E from 
the front profile via the relation 



Q 



long 



hort 



sin{ f ~ 5 arccos B} 
sin{| -f i arccos _B} 



(10) 



with B := |V3{i£;-3/2^ 

once the constant C is known for 

the system. 

The long range shape of the front is determined by 
the smaller root, Qiong, which asymptotically approaches 
the simple law Q — —v/E. Thus the width of an elastic 
front should remain constant during dewetting, in agree- 
ment with recent observations [jisl , but in contrast to 
quasi-Newtonian fronts |p^ , p3| . Since our model predicts 

Qshort to scale as \/~E for sufficiently large we have 
Qiong oc Q^fiort Other parameters in the system are 
kept constant. In fact, for the two curves shown in fig. ^a, 
the ratio of the Qiong is 1.40 ±0.07, while the ratio of the 
Qlhort 1-34 ±0.13. This is clear agreement. 

We thus have put forward a tractable model which 
describes the generic morphologies of liquid dewetting 
fronts observed so far. It correctly accounts for the im- 
pact of the fluid properties upon their occurrence and 
structure. It is thereby capable of extracting information 
on the sub-micron scale rheological properties of the liq- 
uid film from the observed profiles. This may be rendered 



quantitative by determining the numerical constant a, 
which we leave to further work. 
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FIG. 1. Dewetting morphologies observed in various liquid 
polystyrene films beading oflt a silicon substrate. The height 
has been normalized with respect to the film thickness, a) 
Profiles of the rims of holes nucleated in 40 nm thick films. 
The steep profile was obtained with a molecular weight of 
Mw = 101 k, the flatter profile with Mw = 600 k. Each profile 
is superimposed with a double exponential fit as a dotted 
line, the agreement is almost perfect (see text), b) Profile 
obtained with Mw = 2 k. c) Same as b), but with smaller 
film thickness. 

FIG. 2. a) Morphological phase diagram for viscoelastic 
dewetting fronts in the '!;)-plane, according to en. M. b) 
Deviation of the phase boundary according to eq. (|5|) (solid 
curve) from the approximation given by eq. (^) (dashed line). 
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This figure "Figl.png" is available in "png" format from: 
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